UK-04/12 
hep-th/0406029 



On The Problem of Particle Production in c = 1 Matrix Model 



o 
o 

(N 

<: 

(N 

(N 
> 

o 
o 
o 



Oh! 



X 



Partha Mukhopadhyay 

Department of Physics and Astronomy 
University of Kentucky, Lexington, KY-40506, U.S.A. 

El-mail: partha@pa.uky.edu 



Abstract 

We reconsider and analyze in detail the problem of particle production in the time 
dependent background of c = 1 matrix model where the Fermi sea drains away at late 
time. In addition to the moving mirror method, which has already been discussed in 
hep-th/0403169 and hep-th/0403275, we describe yet another method of computing the 
Bogolubov coefficients which gives the same result. We emphasize that these Bogolubov 
coefficients are approximately correct for small value of the deformation parameter. 

We also study the time evolution of the collective field theory stress-tensor with a 
special point-splitting regularization. Our computations go beyond the approximation of 
the previous treatments and are valid at large coordinate distances from the boundary at 
a finite time and up-to a finite coordinate distance from the boundary at late time. In 
this region of validity our regularization produces a certain singular term that is precisely 
canceled by the collective field theory counter term in the present background. The energy 
and momentum densities fall off exponentially at large distance from the boundary to the 
values corresponding to the static background. This clearly shows that the radiated energy 
reaches the asymptotic region signaling the space-time decay. 



Contents 

1 Introduction and Summary 2 

2 Review of Basic Facts 6 

2.1 The World-Sheet Theory 6 

2.2 The Matrix Model 8 

3 The Moving Fermi Sea Backgrounds 10 

4 Bogolubov Analysis 11 

4.1 Bogolubov Coefficients - Scattering Method 13 

4.2 Bogolubov Coefficients - Moving Mirror Method 16 

5 Stress-Tensor Analysis 18 

5.1 Energy and Momentum Densities in Point-Splitting Regularization 18 

5.2 Collective Field Theory Regularization in Static Background 20 

5.3 Collective Field Theory Regularization in Draining Background 21 

6 Conclusion 24 
A Bogolubov Coefficients in Scattering Method 26 
B Geometrical Construction of the In and Out-Modes 27 
C Point-splitting Regularization in the Draining Fermi Sea Background 29 
D Relation between Bogolubov and Stress-Tensor Analysis 31 

1 Introduction and Summary 

Two dimensional bosonic and type OB string theories have non-perturbative dual de- 
scription in terms of the c = 1 matrix model^. Although the CFT description is quite 
complicated [8, 9, 10, 11] the matrix model description is simple. In the singlet sector this 
reduces to the quantum mechanics of an infinite (but fixed) number of non-relativistic 
free fermions in an inverted harmonic oscillator potential. Recently there has been a 
crucial progress in this subject through the understanding of the unstable DO brane [12] 

^The subject is well developed. Various reviews and some original papers can be found in refs.[l, 2, 
3, 4, 5, 6, 7]. 



and its decay on the matrix model side [13, 14, 15, 16]^. Here also the matrix model 
description [16] of this decay turns out to be quite simpler than the BCFT description 
[19]. It is this simplicity of the matrix model description and the fact that it provides 
us with a non-perturbative definition of the theory enable us to probe various issues that 
would have been difficult to address otherwise (see, for example, [17]). Therefore it is 
important to study various other backgrounds, in particular the time dependent back- 
grounds [20, 21, 22, 23, 24, 25, 26, 27], in the matrix model itself to help understanding 
the corresponding world-sheet theories. This might provide us with clues on how to deal 
with time dependent backgrounds in string theory in general. 

Recently in [23], some interesting fermionic configurations in the matrix model were 
viewed as two dimensional cosmological solutions. A particular class of time dependent 
solutions with two parameters were discussed where the envelop of the Fermi sea, having 
the same structure as the static one, has an overall motion on the single particle phase 
space. The motion is such that the Fermi sea fioods in and subsequently drains away 
thus describing creation and the subsequent annihilation of the universe. They have been 
interpreted to be arising from the closed string tachyon condensation [23, 24, 26]. The 
corresponding world-sheet deformation was also suggested. Switching off one of the two 
parameters one reaches two different classes of one parameter solutions where the Fermi 
sea merges with the static one in either the past or the future asymptote. In this case the 
deformation parameter does not really parameterize inequivalent solutions as the change 
in this parameter can be absorbed by time translation^. The solution, hereafter called 
the "draining Fermi sea", in which the Fermi sea starts from a configuration arbitrarily 
close to the static one at early time and drains away at late time has been further studied 
in [24, 25] where the problem of cosmological particle production has been addressed. 
In particular in [25] the Bogolubov coefficients for the particle production have been 
computed explicitly. An approximate analysis of the energy-momentum tensor was also 
performed and it was argued that the contribution to the energy coming from particle 
production and the time evolution of the initial vacuum energy cancel so that the net 
energy remains the same. 

In this paper we study the same questions in more detail. In the static background 

^See [17, 18] for relevant recent works. 

^AU these solutions very much look like the bulk analogs of the boundary rolling tachyon solutions 
discussed in [19]. The one parameter deformations are analogous to the situations where the open string 
tachyon starts from at early time or reaches at late time the top of the potential. In the open string 
context they are interpreted to be either creation or destruction of the D-brane whereas in the closed 
string context they correspond to creation or annihilation of the space-time itself. One crucial difference 
between the two cases is that open string solutions arise from a perturbative instability whereas the closed 
string ones should be considered as different backgrounds altogether [23] because the deformations are 
by non-normalizable modes (there is no perturbative instability). 



the theory is free far away from the tachyon wall. In the draining background the tachyon 
wall moves in such a way that an observer initially sitting in the free region gradually 
enters the strongly coupled region where the notion of the massless particle is lost and 
the theory becomes more complicated. Nevertheless there is a time span during which 
the observer stays in the free region and can observe particle production. We attempt to 
quantify this observation during this particular span of time. Although the deformation 
parameter A corresponding to this background can be changed by time translation we fix 
the origin of time and treat A as an actual parameter. In this time-frame the condition for 
the observer not to enter the strongly coupled region during the whole process of observing 
particle production turns out to be In 1/2A >> 0. We shall see that this condition will be 
essential for computing approximate Bogolubov coefficients for the particle production. 

In the first analysis of the system we study the computation of the Bogolubov coeffi- 
cients in two different methods namely, the "scattering method" and the "moving mirror 
method". The first method is directly related to the wall-scattering'^ in the static back- 
ground and uses crucially the Seiberg bound [35, 36] on the spectrum of primary operators 
in the CFT [8, 9, 10, 11]. The space-time version [3] of this bound is that the in-states 
of this particular scattering problem contains only the right-moving excitations which 
move toward the wall. Similarly the out-states contain only the left moving excitations. 
Application of the same bound in the time dependent case implies that the in and out- 
vacua support only the right and left-moving excitations respectively. In the static case 
the result of the classical wall-scattering is given by a "scattering equation" [33] where 
a right-moving oscillator is expressed in terms of the left-moving oscillators in the form 
of a power series, the leading term being linear and of order g^^. Therefore computing 
the scattering equation up-to the leading order in the time dependent case should give 
us the Bogolubov coefficients. Although this scattering equation is classical, non-trivial 
Bogolubov coefficients, which indicate quantum phenomenon like particle production, can 
be obtained with the following additional input. The motion of the moving Fermi sea on 
the phase-space can be undone by a time dependent canonical transformation which can 
be lifted to a coordinate transformation on the collective field [23, 25]^. Implementing 
this coordinate transformation on the scattering equation in the static background gives 
non-trivial mixing between the positive and negative frequency modes. It turns out that 
this analysis gives approximate Bogolubov coefficients for small A such that In 1/2A >> 0. 

In the moving mirror method, which has already been discussed in [24, 25], we set 
up the whole problem in the framework of the Das-Jevicki collective filed theory [39]. 

^See, for various methods of computing the scattering amplitudes, [28, 29, 30, 31, 32, 33, 21, 34]. 
^Following [25], this can also be viewed as the action of a particular M^oo-transformation [49] on the 
collective field. 



We work at the linearized level of the equation of motion which corresponds to taking 
a formal gs ^ limit^. In this approximation the problem of particle production in 
the moving Fermi sea backgrounds reduces to that of a moving mirror problem with 
reflecting boundary condition. This subject has been studied quite extensively in flat- 
space [42, 43, 44, 45]. In our case, it turns out that the metric is time dependent and 
therefore, as expected, introduces problems in defining the natural modes that should 
correspond to particles. In fiat space examples of moving mirror one uses a certain 
argument of geometrical optics to construct the natural in and out-modes. We generalize 
this construction to the present case at hand. 

The class of moving Fermi sea backgrounds are such that asymptotically the mirror 
approaches the velocity of light. In the draining Fermi sea background, which we con- 
centrate on, this happens at the future asymptote where the mirror moves toward the 
observer'^. In a standard moving mirror example [42] one assumes that at far future the 
mirror approaches a constant velocity less than that of light so that the motion can be 
undone by a Lorentz transformation. In case of a mirror approaching the velocity of light 
toward the observer the whole portion of the future null asymptote is not available for 
defining the modes. In fact the natural out-modes that are constructed following the 
method of geometrical optics form an over-complete set in that region which destroys 
the orthonormality of these modes. It turns out that for the particular mirror-trajectory 
involved in our example the upper bound on the future null infinity is given by In 1/2A. 
Therefore taking lnl/2A >> we make these out-modes approximately orthonormal. 
The Bogolubov coefficients computed using these out-modes are therefore approximately 
correct for small A. As expected from the fact that A acts as a deformation parameter, 
both the above methods give trivial Bogolubov coefficients in the limit A -^ 0. Therefore 
our results are approximately correct for small but nonzero A. 

Next we turn to the analysis of the stress-tensor and go beyond the approximation 
made in [25] (see the last part of the discussion in sec. 6 for more details). Incorporating 
the effect of the non-fiat metric we explicitly compute the vacuum expectation value of 
the stress-tensor by expanding field in terms of the in-modes constructed in the moving 
mirror discussion. Doing the computation in a point-splitting regularization method we 
encounter the usual vacuum ambiguity. Das-Jevicki collective field theory [39] automat- 
ically comes with a particular counter term which fixes this ambiguity (see [25] for a 
recent discussion). This counter term was first obtained in [38] using the collective field 

®The effective coupling, which has gs as an overall factor, increases gradually as we approach the wall. 
By choosing gs arbitrarily small we can increase the region of validity of this approximation. 

^Examples of mirror-trajectories approaching a null line has already been discussed in the literature 
[43, 44] with the exception that mirrors moving away from the observer are usually considered. One of 
the known examples reproduces Hawking radiation [47]. 



method [37] which is formally background independent. It was also showed by Gross and 
Klebanov in their "fermionic string field theory" formulation [40] that this counter term 
is equivalent to the normal ordering directly obtained from the fermionic theory. Consid- 
ering the static case first, we show that in the computation using the mode expansion the 
known result is reproduced by a special type of point-splitting regularization method. We 
generalize this method to the time dependent case and obtain results for both the singular 
and finite parts of the stress-tensor components. We justify this regularization by showing 
that the counter term in the Hamiltonian required to cancel the singular part precisely 
agrees with the collective field theory counter term evaluated at the present background. 
We should mention at this point that our computations are actually done in a metric 
which is simpler than the actual one. This gives rise to a particular region of validity 
only where we achieve the above agreement. The region of validity of our analysis is (1) 
everywhere at early time, (2) large coordinate distances from the boundary at a finite 
time and (3) up-to a finite coordinate distance from the boundary at late time. The finite 
part of the energy density approaches the value corresponding to the static background 
at early time but evolves into something else at late time giving rise to a nonzero radiated 
energy. We show that the energy and momentum densities fall off exponentially to the 
values corresponding to the static background at large distance from the boundary. 

The rest of the paper is organized as follows. We review some of the basic relevant 
points of both the world-sheet theory and the matrix model in sec. 2. The time dependent 
backgrounds discussed in [23] are reviewed in sec. 3. The two methods of computing the 
Bogolubov coefficients are given in subsections 4.1 and 4.2. Sec. 5 contains the stress-tensor 
analysis. Finally, we conclude in sec. 6 where we summarize the basic accomplishments of 
this paper and compare our stress-tensor analysis with those preexisted [24, 25, 26] in the 
literature. The appendices contain some of the technical details. We discuss the relation 
between the Bogolubov and the stress-tensor analysis in appendix D. 

2 Review of Basic Facts 

Here we touch upon the basic features of the world-sheet and the matrix model descrip- 
tions relevant for our discussion and spell out the dictionary of the duality. 

2.1 The World- Sheet Theory 

In addition to the usual time-like scalar X^, the matter part of the world-sheet theory 
contains only one space-like scalar X^ whose action is given by a particular limit of the 



Liouville action [8, 9, 10, 11], 

Sl = ^ J d^z [dX'dXi + 271 fioe^"'"') , (2.1) 

with a background charge Q = b + 1/b and central charge c^ = 1 + 6Q^. The space of 
normahzable states is given by the collection of the conformal families corresponding to 
a one parameter family of primaries given by* , 

Vq+,^ = e(«+™)^' , 0<w<oo, (2.2) 

with conformal dimension (Q^ + w^)/4. Notice that w is restricted to be positive, the so- 
called Seiberg bound [35, 36]. In fact the operators with negative values of w are related 
to the above ones by the following reflection equation, 

Vg+iw = R{u])VQ-iw , 
o.,„^ ^ ,^;, Til + tw/b)Til+twb) 

^' ^ T{l-iw/b)T{l-twb)' 

/i = 7r/io7(&') , l{x) = f^^ • (2-3) 

The world-sheet theory S'xi relevant for string theory is obtained (at a' = 1) by taking 
the following limit, 

6^1, ;Uo -^ oo , /i fixed . (2.4) 

It turns out that it is only the combination fi/gs which appears as a parameter in all 
the amplitudes so that there is only one parameter of the string theory which we take 
to be Qg. BRST analysis of the world-sheet theory shows that there is only one space- 
time field-theoretic degree of freedom, namely the "massless tachyon" given by the vertex 
operators, 

V^ = cce-™^V2±,^ = cce-^^-^^'^^^'^e^^' , w>0. (2.5) 

The space-time interpretation of the reflection equation (2.3) or the Seiberg bound w > 
in the above equation goes in the following way (see, for example [3]). The effective 
coupling Qeff = QsG^^ falls off to zero at large negative x^ where we have a free massless 
particle [36]. Because of the Liouville interaction- wall in (2.1) any right-moving pulse 
created in the free region necessarily results into left-moving pulses due to the scattering 
from the wall. Therefore in this scattering problem the in and out-states always contain 

*Note that this form of the vertex operator is vahd only at large negative X^ where the interaction 
term in (2.1) is negligible. 



only the right (V^) and left-movers {V~) respectively. In the second quantized theory 
these states are created by harmonic oscillators a™''' and a'^*^ respectively. We take the 
following normalization. 

K,€^] = [aZ\ay']=S{w-w'). (2.6) 

Then one concludes from (2.3) with the proper limit (2.4) that to the leading order in 
Qs, a^ is proportional to a°"* where the proportionality constant is simply a w dependent 
phase. 

2.2 The Matrix Model 

In the singlet sector, the c=l matrix model reduces to the quantum mechanics of an infinite 
(but fixed) number of non-relativistic free fermions in an inverted harmonic oscillator 
potential given by — x^/2. Qs plays the role of the Planck's constant h of the fermion 
theory. The string theory Hamiltonian is given by the second quantized Hamiltonian 
of these fermions with an additional factor of I/Qs [33]. The closed string background 
discussed in the previous section corresponds to a particular classical state where all the 
single particle phase-space trajectories below energy —1/2 are filled by fermions. For the 
bosonic case, which is under consideration in this paper, only those trajectories which are 
on the left side of the potential are filled so that the envelop of the (static) Fermi sea is 
given by, 

{x + p){x-p) = 1, x<-l. (2.7) 

Let p±{x,t) be the values of p at the fiuctuating upper and lower edges (fig.l). Then the 
fiuctuations ri±{x,t) are defined as, 

P±{x, t) = po±{x) + 77±(x, t) , (2.8) 

where the background is given by, 

Po±{x) = ±Po{x) , Po{x) = Vx^-l , x<-l . (2.9) 

Using Polchinski's bosonization in [33] we define the fiuctuations of the scalar field and 
the conjugate momentum, denoted by (j){x,t) and 7i^{x,t) respectively, 

r]±{x,t) = y/ngs -7if{x,t) ± dx(t>{x,t) . (2.10) 

X = —1 acts as a boundary where satisfies a Dirichlet boundary condition^: (f){x = 
— l,t) = 0. The static Fermi sea background (2.9) is given, in terms of the collective 
^This can be derived from the constraint that the total number of fermions be fixed. 
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Figure 1: The "static Fermi sea" background which is dual to the closed string background 
with linear dilaton and static tachyon wall. The fluctuation of the upper and lower edges 
(which propagate along the direction of the arrows) give rise to the massless tachyon field. 



variables, as, 



(9^0f)"^^(x,t) 






tTtcs) (x, t) = . 



(2.11) 



Defining the new coordinate q, 



X 



cosh q 



-oo < g < , 



(2.12) 



so that Pq = I sinhgl, Podx = &„, the linearized equation of motion takes the form. 



id^ - dl 



0, 



(2.13) 



where 4>{q,t) = (f){x,t). It turns out that at g << the non-linear parts of the equation 
of motion reduce to an infinite series of interaction terms with Qeff = QsG^^ (g = g + ln2) 
as an effective coupling. In this region the fluctuation 0(g, t) can be mode expanded as. 



0m (m) + (j^outiv) , 

"^ dw 



[u 



^out{V 



y^TTW 

°° dw 



in —iwu 



^/4' 



TTW 



out —iwv 



+ h.c. 
+ h.c. 



(2.14) 



where u = t — q and v = t + q are the usual null coordinates. The oscillators a^ and 
a^"* have the same normalization as in (2.6). The scalar field S{q,t) introduced in [33], 



is related to (p{q,t) at q « in the following way: S{q,t) = (f){q,t). Therefore the 
oscillators, say a^ and a°"*, of S{q,t) are related to that of (j){q,t) as follows: a^ = 
gj«)in2^m^ ^out _ ^-iw\n2^out ^ Using the result obtained in [33], we get the following 
scattering equation, 

a- = e^-i-4ar + 0((?.). (2.15) 

The oscillators of the scalar field directly coming from string theory are related to the 
above oscillators through the leg-pole factors [28, 41, 32] with the additional phase factor 
described above. 

-in _ -*- V-"^) -iwln2 in -out _ -*- \ "^"^J iw hi 2 out /cy -r (-\ 

The above derivation of the collective field theory is classical. Therefore in the process 
of quantization one encounters the usual vacuum ambiguity which has to be fixed by 
a normal ordering prescription. In the string field theory of Gross and Klebanov [40], 
which was shown to be equivalent to that of Das and Jevicki [39], this automatically 
came from the normal ordering of the fermionic theory. In Das- Jevicki field theory this is 
incorporated by a particular counter term first computed in [38] using the collective field 
method [37]. In our notation this counter term is given by, 

AH = — ^-= / dx dx(j){x, t)dxdxi In \x — x'\x=x' ■ (2-17) 

2y7r J-oo 

According to the collective field derivation [38] this counter term is background indepen- 
dent ^°. 



3 The Moving Fermi Sea Backgrounds 

We shall now consider a particular class of two-parameter time-dependent closed string 
backgrounds discussed in [23]. On the matrix model side this is described by a moving 
Fermi sea whose envelop is given by, 

(x-hp + 2A+e*)(x-p-h2A_e"*) = 1, x < x^{t) , ^^(t) = -(1 -h A+e* + A_e"*) . 

(3.18) 

The envelop of the Fermi surface is an hyperbola which moves like a rigid body following 
the hyperbolic trajectory of its centre: a;o(t) = — (A+e*-|-A_e~*) , poit) = — (A+e* — A_e^*). 

-'^"l am thankful to S. R. Das for discussion on this point. 
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Just like as described in the previous section the fluctuation of this Fermi surface also 
gives rise to a (1 + l)-diniensional scalar field 0(t, q) which, at any given instant of time, 
behaves like a free massless field at a sufficiently large negative q. On the phase-space 
this configuration can be transformed into the static sea by the following time-dependent 
canonical transformation^^, 

X = a; + A+e* + A_e"* , P = j) + A+e* - A_e"* , -oo < X < -1 . (3.19) 

Following the usual method as described in subsection 2.2 a (1 + 1) dimensional scalar 
field theory can be constructed from the static configuration in the XP-plane. We shall, 
in this case, use upper case symbols for various objects discussed in the previous section. 
We relate the space-time of this field theory with the original (matrix model) space-time 
(x, t) by the following coordinate transformation, 

X = a; + A+e* + A„e-* , T = t. (3.20) 

Following (2.12) we define new coordinates Q and q as, 

X = - cosh Q , -oo < Q < , 
X = — coshg , — oo < q < qm(t) , 
q^it) = -cosh-^(l + A+e* + A_e-*) , (3.21) 

such that, 

cosh Q = cosh q - (A+e* + A_e"*) . (3.22) 

At large negative Q this reduces to Q = q — ln(l — 2A+e^ — 2A_e^") such that, 

U = u + ln(l - 2A+e^ - 2A_e~") , V = v- ln(l - 2A+e^ - 2A_e-") , (3.23) 

where as usual we have defined: U = T — Q, V = T + Q. The scattering equation for the 
oscillators A^ and A°"* of the scalar field ^{Q,t) will then be given by the same equation 
as (2.15), 

4 Bogolubov Analysis 

For our explicit analysis we shall consider the draining Fermi sea background in which 
case we have in mind all the equations of the previous section with A+ = A, A_ = 0. 
On the matrix model side the wall scattering of the massless particles can be computed 
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Figure 2: In the draining Fermi sea background the observer initially sitting far away 
from the wall at go << gradually enters the strongly coupled region. 



[33] very easily by exploiting the fact that the classical trajectories of the fermions on the 
phase-space are exactly known. Given an incoming pulse created on the upper edge of the 
Fermi sea at g << one can compute the deformed outgoing pulse that re-emerges, after 
some time on the lower edge without worrying about the complicated self interactions 
that the pulse goes through. In the draining Fermi sea background, as time evolves the 
wall moves toward left (fig. 2) so that the observer gradually enters the strongly coupled 
region. We want to study the system over a time span during which the observer remains 
in the free region. If a pulse is sent in by the observer at a time t_(<< 0) which is taken 
to be of the same order of the observer's position qo{<< 0) then the time at which the 
pulse is received back is given by, 

t+ = -q^ - ln(l + 2A) . (4.25) 

The boundary of the space qm{t), which is zero at t = t_, is given, at t = t+, by, 

2A \ 



Hm\^- 



ln2A-t+ + 0(l/Ae*^ 



go -In 



1 + 2A 



(4.26) 



"'^^As discussed in [25] this can also be viewed as a particular Woo transformation. 
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where we have taken e'^° « y^2A- To keep Qo in the free region at t = t+ we should have 
Qo — (lm(t+) « 0. These two conditions can be written as, 

-go » In 1/2A » . (4.27) 

We shall see that to compute the Bogolubov coefficients we need to work with sufficiently 
small A so as to be consistent with the above equation. Then given a small value of A it 
also gives an estimation for the position of our observer. 

4.1 Bogolubov Coefficients - Scattering Method 

Let us now proceed to compute the Bogolubov coefficients. These are the coefficients of 
the linear expansion of the out-oscillators in terms of the in-oscillators where a mixing 
between the positive and negative frequency modes occur. One may, therefore, wonder 
if the leading order scattering equation obtained by the method of [33] can give these 
coefficients. In spite of the fact that the analysis of [33] gives only the tree-level results, 
with an additional input of coordinate transformation it is indeed able to capture particle 
production which is a quantum effect. At late time, the natural oscillators in the draining 
Fermi sea background are a^"*, whereas at early time these are given by, 

a™ = At" = e™^°^A^"* + 0{gs) (4.28) 

where the ffist equality is due to the fact that at early time the draining background 
merges with the static one and in the second step we have used the scattering equation 
(3.24) in the static background. The Bogolubov coefficients, therefore, effectively relate 
a°"* and A°"* oscillators. 

a^^' = / dw' \a{w, w')Al^' + f3{w, w')AZ'^ • (4.29) 





These coefficients can then in turn be computed by using the transformation, 

(9^0 = (9^f/(9t/$ + d^,Vdv^ , (4.30) 

and setting du^ = at late time which is the main use of Seiberg bound in the present 
analysis. For small A we get the following results (see appendix A for details), 

nnr /•lnl/2A rl-r , , , 

V W J-oo 271 

f3{w,w') = /^/'°'^'^^e™^+™'^-(^) , (4.31) 

V w J-oo zn 
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Mirror trajectory 
t 



Figure 3: The asymptotic null infinities are given by, X^ : (—00 < u < 00, v = —00), 



II 



u 



00, — 00 < t> < f„ 



1/ ln2A) and the mirror trajectory is given by, u = f\{v). 



The mirror approaches the null line v = Vm = lnl/2A at late time. Taking A ^ pushes 
Vm to infinity resulting in a trivial trajectory. 



where the function f\{x) is given by, 

/^^(a;) = a; - ln(l - 2Ae^) , -00 < a; < In 1/2A . 
Several comments are in order. 



(4.32) 



1. The results (4.31) resemble the standard expressions for the Bogolubov coefficients 
in a moving mirror problem in fiat-space with mirror-trajectory (fig. 3), 



u = fx{v) . 



(4.33) 



Notice that this is not the actual trajectory given by qm{t). We shall understand 
the physical origin of the function fx{x) in the above equation when we analyze the 
moving mirror method in the next section. 



2. Notice that, 



hm a(w, w') = 6(w — w') , lim 3(w, w') = . 



(4.34) 



indicating no mixing of modes, hence no particle production. This corresponds to 
the static background as the trajectory trivializes in the above limit. 
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3. We should emphasize that the above Bogolubov coefficients are approximately cor- 
rect for small but non-zero A. The small A limit is necessary for the following 
approximation which has been applied in our derivation (appendix A) to invert a 
certain equation. 

/•lnl/2A rl-r 
J-oo ZTC 

~ 6{w) , for small A . (4.35) 

4. Explicit evaluation of the integrals in (4.31) give, 



a 



2yw' zn \ w T(iw) 



2\ w' 2-n \ w T[iw) 

(4.36) 

• Keeping in mind that the frequencies are always positive so that we can set 
5{w + w') = one recovers the trivial result (4.34) in the limit A ^ 0. 

• For small but non-zero A the second term in each of the above expressions 
are highly fluctuating. The non-trivial feature that the A dependence appears 
only in a phase makes it possible that certain physical quantities in which 
this phase is canceled have A independent results (we shall see an example 
below). This is consistent with the fact that A can actually be changed by 
time translation and therefore a typical time-averaged quantity should not 
be sensitive to this parameter. Local quantities which have to be given at a 
particular time do depend on A. Although it is not clear how to make sense of 
these fluctuating objects in that case we argue in appendix D that the energy- 
momentum analysis give smooth A dependent functions which can be related 
to some integrals of these Bogolubov coefficients. 

The above Bogolubov coefficients relate the oscillators a^"* and A^"*. We may relate the 
"string theory" oscillators a^"* and a™ by introducing another set of coefficients a^{w, w') 
and P^{w,w'), 



a°"* = / dw' 



/O 



a^{w, ti;')C + P'{'^. ^')Cl • (4.37) 



Since at early time the present background approaches the static one eq.(2.16) may be 
used to relate the "string theory" and "matrix model" oscillators. But it is not clear if 
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eq.(2.16) is still valid at late time. Assuming that this is the case one can derive the 
following expressions, 

a^{w,w') = V{w,w')a{w,w') , 

(3^lw,w') = V{w,-w')(3{w,w') , (4.38) 

where V{w,w') is simply a pure phase, 

V{w,w') = e-(— ')^"^ ^^p"^y"\ (4.39) 

so that the result for the total number of out-going particles between the frequency range 
w and w + dw remains the same, 

N{w) = / dw' p^{w,w') = dw'\P{w,w')f , 





oo 





dw' sinhTTW ,, ,„, 

. (4.40) 

47r(w + w') sinh7rw'sinh7r(w + w') 

This is an example of a A independent physical quantity mentioned above. None of our 
expressions are infrared regularized. Apart from this difficulty the above integral is well 
behaved. In fact, with an infrared cutoff in the integral, N{w) is ffnite for all w. 

4.2 Bogolubov Coefficients - Moving Mirror Method 

We have seen in subsection 2.2 that at the linearized level the static background corre- 
sponds to a massless scalar on a half-line with a ffat metric. In this approximation the 
problem of particle production in the draining Fermi sea background reduces to that of a 
moving mirror problem in a non-ffat metric [24, 25] . The metric Qfj^^, in (t, g)-coordinates 
can be obtained by transforming the metric, G^i, = diag(l,-l), which corresponds to 
(T, (5)-coordinates under the coordinate transformation (3.22) with A+ = A and A_ = 0, 

_ 1 / sinh^ g — 2Ae* cosh g Ae*sinhg\ {AA^'\ 

^^" ~ (coshg-Ae*)'-l V Ae*sinhg - sinh^ g J ^- > 

At large time this reduces to, 

1 ( 1 - 4Ae^ -2Ae^ \ 

^^- - (1 _ 2Ae-)2 - 4e29 1^ -2Ae^ -1 ) ^^'^^^ 

Notice that at late time eqs.(4.26, 4.27) give: qm(t+) ~ go + In 1/2A << 0. Therefore for 
any g < qm{t) at which (1 — 2Ae^) is ffnite, e^^ is negligibly small. In this approximation 
the metric turns out to be, 

1 / 1 - 4Ae'' -2Ae^ \ 

^^- - (1 _ 2Ae-)2 l^ -2Ae^ -1 ) ' ^^'^'^^ 
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Notice that even at finite time the metric (4.41) reduces to the above one at large |g|. This 
form precisely corresponds to the coordinate transformation (3.23) with A_|_ = A, A_ = 0. 
At early time when v -^ — oo, g^j^^ -^ G^^. This implies that the metric (4.43) which looks 
simpler than (4.41) captures the correct coordinate systems in the following region of 
validity: (1) everywhere at early time, (2) large coordinate distances from the boundary 
at a finite time and (3) coordinate distances from the boundary such that (1 — 2Ae*') is 
finite at late time. To avoid complications we shall replace the metric (4.41) by (4.43) for 
all space-time in our analysis. In this simplified model the mirror-trajectory is given by, 



u = p(v) , piy) 



21n(l-2Ae" 



(4.44) 



Now the question is: given the metric (4.43) and the mirror-trajectory (4.44) what are the 
natural in and out-modes? In appendix B we have generalized the argument of geometrical 
optics for constructing these modes from flat-space to the present case. The results are. 



tout 



1 



V4 



TTW 



Q-iwS, _|_ g Q-iwfx{v) 



^-iwgxiO _|_ g ^-iwv 



(t) - g) , w>0 

d{qra{t) - g) , W>f) 



(4.45) 



where ^(-u, f ), given in eq.(B.6), is one of the null coordinates in the present metric which 
replaces u in the flat case. The other null coordinate v remains the same as in the fiat case. 
The functions f\{x) and gx{x), which are inverse of each other f\{g\{x)) = x, are given in 
eqs.(4.32) and (A. 2) respectively. Notice that, as discussed in appendix B (eq.(B.8) with 
C = V and f{x) = fx{x)), the function f\{x) characterizes the trajectory in terms of the 
null coordinates (^,f), which is the key reason for it to appear in the expressions (4.31) 
for the Bogolubov coefficients. 

There remains the question of orthonormality of these modes, a property which is 
needed to compute the Bogolubov coefficients. We shall now discuss this issue. The 
"time" of our problem is inherited from the matrix model to be the coordinate t. It 
can be checked that this is a reasonable candidate for a "time function" in the present 
case. Therefore the corresponding space-like surface S is given by the constant t surface 
coordinatized by q. Then the inner-product is given by^^ 



rqmit) 



n,< 



—I 



J —oo 

1 + 2Ae'' 



oo 

qm(t) 






dq 



oo 



1 - 2Ae^ 



01-9702 + 01'9r02 



(4.46) 



12 



Since g — dety^i/ — gii, the lapse function N [46] is simply 1 in our case. 
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where we used the usual notation: Xd*^Y = XdY — {dX)Y. Using this inner product one 
can exphcitly check that the in-modes (4.45) are always orthonormal, but the out-modes 
are orthonormal only for small A when eq.(4.35) is valid. 

(CO = '^(^1 - ^2) , {K^<I>Z*) = , always , 

(C/, C*) = '^(^i - ^2) , (0:^/, 0:f ) = , for small A . (4.47) 

In the computations of (0°"*,0^"*) and (0^"*,0^"**) the function I\{w) (see eq.(4.35)), 
with proper arguments, appears at various places. In the small A limit, which is relevant 
for our discussion (see eq.(4.27)), one uses eq.(4.35) to achieve orthonormality^^. Recall 
that the computation of the Bogolubov coefficients in the scattering method, described 
in appendix A, also went through only for small A. 

In the present method one can compute the Bogolubov coefficients by computing 
certain inner-products of the in and out-modes. Defining 

a{w, w') ^ {<pz, 0r ) , f3{w, w') ^ {<i^:r, 0r) , (4.48) 

we get the same results as reported in eqs.(4.31). 

5 Stress- Tensor Analysis 

So far we have studied the problem through the Bogolubov method. By construction this 
method addresses the question of how many particles are created in a time dependent 
background. A field theory also comes with a natural definition of an energy-momentum 
tensor. In a time dependent background one can ask how the vacuum expectation value 
of the energy and momentum densities evolve with time. This is the question that we 
shall be interested in. 

5.1 Energy and Momentum Densities in Point-Splitting Regu- 
lar izat ion 

The unregularized energy-momentum tensor is given by, 

1 
T^r, = (9^0(9^0 - -gt^uT , 



^^The key reason behind the lack of orthonormahty for a finite A is that T^ is not a full line, rather it is 
truncated due to the presence of the horizon at w = Wm = In 1/2 A (fig. 5). Notice that for the cases where 
the mirror reaches a constant velocity less than that of light the whole portion of T^ is available and 
f/*""* constructed in appendix B are orthonormal. This problem also exists in the usual moving-mirror 
analysis in flat space [44]. 
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T 



g'^d. 



(5.49) 



The canonical Hamiltonian and the total momentum are defined to be^^, 



H 



V„ 



gm{t) 



dqV^ Vo , P 



qm{t) 



dq^/—g V 



1 ) 



2 r^T^, 



T, 



flU) 



(5.50) 



Working in the Heisenberg picture we always take the vacuum to be the in-vacuum, 
the one that is annihilated by a^'s corresponding to the modes 0^ given in eqs.(4.45). 
The vacuum expectation values of the above operators can simply be obtained by mode 
expanding the scalar field. We use the in-modes and adopt a symmetrical point-splitting 
regularization to separate out the singular and the regular parts. 



ai 



Sa 

y ' 



02 



Sa 



(5.51) 



The details of this procedure are given in appendix C. Below we quote the final results 
for the energy and momentum densities. The region of validity of these results is same as 
that mentioned below eq.(4.43). 



S 



-9 (Vo) 
1 



4:7r{6v) 



= £s 
1 + - 



£ 



R 



hi 



(7 - /^a)^ 



^{h^){l-hl) 
27r(7 - h^f 



+ 



{h',f 



167r(7 — hxY-" 



S{fx,v) 
247r 



(5.52) 



V 
Vs 



-9 (Vi) 
1 



^{hx){l + hx '' 
27r(7 - hxY 




S{fx,v) 

247r 



(5.53) 



where the subscripts S and F refer to the singular and the finite parts. The point-splitting 
is given by 6v and the limit 6v —^ is understood. The path along which the points are 
split is given by the function 7(-u, f). 



6u = 'y{u,v)6v 



(5.54) 



^^The synimetrization between the indices /i and v in the third equation docs not have any effect on 
the finite part of the vacuum expectation value. Certain hnearly divergent terms (C.19) obtained in the 
point-sphtting regularization get canceled because of this symmetrization. 
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Using the metric (4.43) one can check that the path is null whenever 7 = h\^ which is also 
expected from the results (5.52, 5.53) [45]. The functions hx{v), ip{h\) and the Schwarzian 
S'(/, x) are given by, 

hx{v) = ~d,i = 2\e^rM, 

,/, N hx h\ h\ 

ih(hx) = — + — + —, 
^^ ' 24 8 12 ' 

Off ^ n^) 3/r(a;)y 

where the primes denote derivatives with respect to the argument. Notice that since we 
have used the in-modes (which are always orthonormal), formally the results (5.52, 5.53) 
do not need a small A approximation. But as explained at the beginning of sec. 4 the 
linearized approximation to the equation of motion may not be good at a finite A. 

5.2 Collective Field Theory Regularization in Static Background 

Physical conclusion can not be drawn from the results (5.52,5.53) as we have not yet 
fixed the subtraction scheme. Since we simply have a scalar field in a background metric, 
normally this procedure will be ambiguous. But as we mentioned toward the end of 
subsection 2.2 that being a dual description of a string theory, the collective field theory 
Hamiltonian automatically comes with a particular counter term (2.17) which fixes the 
subtraction scheme. Below we shall find a way to incorporate this subtraction scheme in 
the present method of computation through mode expansion. This is easier to do in the 
static background which we discuss first turning to the draining background in the next 
subsection. 

Since we are dealing with a free {qs -^ 0) theory the renormalization is only at the 
one loop level. Therefore we simply replace dx(j){x,t) in the counter term (2.17) by the 
background solution (2.11). 

AH = — / dxPQ{x)dxdx' In \x — x'\x=x' , (5.56) 

27r J-00 

where the function P{){x) is defined in eq.(2.9). Moreover, the vacuum expectation value 
of the renormalized energy density (in q space) is given by [39, 40]^^, 

^r = -^S{q,x), (5.57) 



15, 



SCoU _ _i/247r 



At large negative q, where the effective coupling is small and the present analysis is reliable, one gets, 

>i 

R 
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where the function q{x) is given in eq.(2.12). At early time one gets from (5.52), 
Now consider the following point-splitting regularization, 

Noticing that Po{xY = l/{q'{x)Y we write, 

^ - -"^^^"^ 2n{q{x,) - q{x2)y ' ^'-'^^ 

Then using the following identity, 

2 = /r^\^ + g'^l-^?^) + 0[dx) , (5.61) 



(ir(a; - fe/2) - ir(x + (5a;/2))' {dxf 6 

for a function K{x) we arrive at. 

The first term gets canceled precisely by the counter term in the collective field theory 
Hamiltonian (5.56) so that the renormalized energy density is same as that given in 
eq.(5.57). 

5.3 Collective Field Theory Regularization in Draining Back- 
ground 

Here we generalize the regularization (5.59) to the draining Fermi sea background. Given 
the solution (2.11) for the static background we may write the regularization (5.59) in 
terms of this solution which can then be naturally generalized to, 

f(9^0l^Vx,t)f 

cColl ^ cColl _|_ cCoU ^ V ^^0 ^ ' ' J c 

''°^''^' Ss + Sn, (5.63) 



Po(Xi)Po(X2 

where we have used the solution for the draining Fermi sea, 

P^{X) 



dj^\x^t) = ^^ (5.64) 
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which can be obtained following the discussion in section 3 and noticing that d.^ 

?c 

'5 



To compute the singular and the finite parts £g°^^ and £p°^^ respectively one uses the 



following results^^, 
Po{Xf 



Po{X,)Po{X2) {6vy 



Poixy 



(^ + >,X) 



5v 



6v 



i(^ 



7 



2h,f 



^ + 2hx)i'ihx)i5vy + Oi5v' 



(5.65) 



Notice that Xi and X2 are not symmetrically separated from X. Therefore using the 
identity (5.61) directly gives the Schwarzian term evaluated at X which sits midway 
between Xi and X2. But since this term is not multiplied by any singular term it is 
finally evaluated at X in the coincidence limit. The results for Sg"'''' and Sp"'''' turn out to 
be, 
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2h,? 
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hi 



^(/ia)(7'-27/^a + 1) 
7r(7-/iA)2(l-7 + 2/iA 



Sf 



il-hx)\ 
A similar analysis for the momentum density gives the following results, 



(5.66) 
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(5.67) 



The Schwarzian terms in the above equations can be obtained by using the result, 



P^{X)SiQ,X) 



+ 



2 ■ 2P2(X) ' 

and noticing that in the region of validity of our analysis so that we have, 

77(7, /ia) ij{hx){-f^-2-fhx + l) 



(5.68) 
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247r 7r(7 - hx^il - -f + 2hx 
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R 



-'^^Various coordinates like X, Q can be recalled from eqs.(3.20, 3.21, 3.22, 3.23) with A+ = A, A_ = 0. 
As argued below eq.(4.43) we can take eq.(3.23) to be correct for our simplified model. This gives the 
second equation. 
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Mirror trajectory 
t 



Figure 4: y = \v — lnl/2A| gives the coordinate distance between the projection of an 
event on X^ and the boundary on X^. At late time, since qm{t) ~ — t + In 1/2 A, y is same 
as the coordinate distance from the boundary. 
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(5.69) 



Considering a spatial splitting {'~f{u,v) = —1, which we shall justify below (5.73)) we see 
that both the above results depend only on v. In fact defining, 



y 



-(t;-lnl/2A) , 



hiv) = hx{v{y)) 



l-e-y 



(5.70) 



the results in (5.69) can be expressed entirely in terms of y (fig.4). The results turn out 
to be quite simple for 7('U, w) = —1, 
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247r 487r 
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(5.71) 



so that at large y both the above quantities fall off exponentially to the values corre- 
sponding to the static background. 

Notice that the final results (5.71) very much depend on what regularization we work 
with. We shall now justify our regularization (5.63) by showing that in the region of 
validity of our analysis, it is consistent with the collective field theory counter term. The 
original derivation [38] of the collective field theory counter term (2.17) was formally 
background independent. This would imply that the counter term for the draining Fermi 
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sea background will be given by, 

1 fXm(t) 

= 7—/ dx Po{x + \e) dxdx'ln\x -x\^=cc' . (5.72) 

ZTT J-oo 

In the regularization (5.63) we have an undetermined function '~f{u, v). As long as we take 
lim^^_oo7(M,'y) -^ —1 we recover both the singular and the finite parts of eq.(5.62) from 
eqs.(5.66). But for a Hamiltonian formulation it will be more natural to consider a spatial 
splitting. This requires us to take for all u and v, 

-f{u,v) = -l. (5.73) 

With this one can compute Sg"'''' given in eqs.(5.66). To cancel this particular value of 
gCoii Qj^g needs the following counter term, 

{AH,U = ^J_ dx I p^(^) I -j^ ■ (5.74) 

It is straightforward to check that the factor in the curly bracket simply reduces to 1 in 
the region of validity giving the same result as (5.72). 

6 Conclusion 

The problem of particle production in the draining Fermi sea background of the two 
dimensional string theory was considered in [25] (see also [23]). In this paper the same 
problem has been reconsidered and analyzed in more detail. The basic accomplishments 
are the following: 

1. Physical arguments have been given to establish that the Bogolubov analysis is 
sensible only in a limit where the deformation parameter A is small but nonzero. 

2. An additional method of computing Bogolubov coefficients, namely the scattering 
method (which is very reminiscent of Polchinski's wall scattering method [33] in the 
static background), has been presented. 

3. The moving mirror method of computing the Bogolubov coefficients has been elab- 
orated through the explicit construction of the natural out-modes. We show that 
the small A limit that has been considered on physical grounds is required to hold 
in order for both the above methods (which are quite independent) to work. 
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4. Time evolution of the energy-momentum tensor in the hnearized approximation 
of the equation of motion has also been studied in detail. In [25] the coordinate 
transformation (3.23) (with A_|_ = A , A_ = 0) was taken to be, 

U = u , V = v- ln(l - 2Ae^) , (6.75) 

which is a conformal transformation. This, in fact, seems to be a good approxima- 
tion-'^'^ on X^^ (where u -^ oo) at a finite value of ln(l — 2Ae^). The present work 
attempts to go beyond this approximation by incorporating the effect of the second 
term in the expression for U in eq.(3.23). This term takes the transformation off 
the conformality giving a non-trivial effect on the metric (see eq.(4.41)). 

The present computations have been done with a metric (4.43) which is simpler than 
the actual one. This procedure gives rise to a reduced region of validity given by, 
(1) everywhere at early time, (2) large |g| at a finite time and (3) finite (1 — 2Ae^) at 
late time. Our computation is done in a particular regular izat ion which leads to the 
certain singular term that is required to cancel the collective field theory counter 
term for the present background. As expected, this cancellation takes place only in 
the region of validity mentioned above. 

The final results show that the energy and momentum densities have non-trivial time 
evolution. They start up with constant values at early time but grows exponentially 
at late time as one approaches the boundary from a large distance. This clearly 
indicates radiation of energy due to the space-time decay. This is different from the 
result found in [25]^^ where the densities remain constant as time evolves. 

Similar remarks apply to the result found in [24, 26]. The only difference between the 
treatments of [24] and [25] is that in [24] the initial value of T^^ (which corresponds 
to the static background) was taken to be zero as opposed to the case of [25]. As 
a result [24] also finds exponential fall-off of the energy density. Although this is 
similar to the result found in this paper, the reasons are different. 
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A Bogolubov Coefficients in Scattering Method 

Here we give the details of the computation of the Bogolubov coefficients deffned in 
eq.(4.29). Using eq.(4.30) with du^ = 0, the coordinate transformation (3.23) with 
A+ = A, A_ = and the mode expansion for the field we get, 

-iwgxiV) 



dwJw 



J^out^-iu,V _ ^_^_ 



dwJw 



,out 



1 + 2Ae^ 



— h.c. 



(A.i; 



where we have introduced the function, 

^^(x) = a; - ln(l + 2Ae'') , -oo < a; < oo . (A. 2) 

Using the expansion (4.29) in eq.(A.l) and performing an indefinite integral over V on 
both sides one gets, 

J_^-^n>v ^ f^ dw^ r (^,^ ti;)e"^"''s^(^) + f3* {w' , w)e'^'^'^''^] (A.3) 

^/w Jo \/w' L ^ 

Introducing the function fx{x), the inverse of g\{x), as in eq.(4.32) we write the identity 
(A.3) as. 
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Multiplying both sides by e*™"^ and integrating over the valid range of x (see eq.(4.32)) 
one gets. 
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(A.6) 



-oo 27r 

which, in the small A limit, approaches 6{w). Approximating Ix{w) by S{w) we get the 
coefficient a{w,w) from eq.(A.5). Similarly multiplying both sides of eq.(A.4) by e"™^ 
and following the same steps one evaluates the coefficient P{w,w). The results that one 
gets are given in eqs.(4.31). 



26 



B Geometrical Construction of the In and Out-Modes 

In a moving mirror problem in flat space the natural in and out-modes have the following 
geometrical optics interpretation. Each of the two kinds of modes is a superposition of an 
incident and a reflected waves. The incident wave in an in-mode takes the simple plane- 
wave structure on X^, but the reflected one is complicated on X^. For an out-mode the 
incident part is complicated on X^ so that after reflection it takes the simple plane-wave 
form on X^. This means that if the metric (4.43) were flat the trajectory (4.44) would 
have corresponded to the following natural modes, 



I out 
w 



r-, . ^^ ^(gm(t)-g) , w>Q, 



where p~^{x) is the inverse function of p{x) and s = =F correspond to the Dirichlet 
and Neumann boundary conditions respectively. We now try to generalize the above 
construction to our case through the following steps, 

1. Find the two null directions R±{u,v) which can be found, up-to overall factors, by 
solving the following equations, 

RUu,v)R±^{u,v) = 0. (B.2) 

We get, 

R^r^il , 1 - 4Ae^) , i?^ ~ (1 , - 1) . (B.3) 

These give the structure of the light-cones (fig.5) in the given space. Since the 
metric becomes flat as f ^ — oo, the light-cones are stretched to 90° cone-angle 
in that region. They shrink as v increases and become completely squeezed at 
V = Vm = lii 1/2 A where a horizon is formed^^. 

2. Find the null rays 

C{u,v) = const. C{u,v) = const. , (B.4) 

i.e. the the curves which have R'^{u,v) as tangent vectors respectively. 

R^d.^u, v) = 0, R'td.Ciu, v)=0. (B.5) 



^^Notice that this is true only for the simphfied metric in (4.43) and not the actual one given in 
eq.(4.41). For the actual metric the light-cones are completely squeezed at the trajectory of the mirror 
itself where the metric is singular. 

27 




Figure 5: Curves with dot-style 1 and 2 are the null-rays whose tangent vectors are given 
by i?+ and R^ respectively. Light-cones are same as that of a flat space on X^. They 
shrink as v increases and are completely squeezed at w = f^ = In 1/2A where R'^ = R^. 



The two null coordinates are found to be, 



^(M,t;) = M + ln(l-2Ae^) , ({u,v) = v 



(B.6) 



3. Generalize eqs.(B.l) by the following, 

1 



,in 
w 



I out 



1 



V4 



TIW 



e ^ + s e 



-iwfiO' 



+ s e 



'iwC, 



d{qm{t) -q) , w>0 , 
(t) - g) , w>0, 



(B.7) 



where the functions f{x) and g{x) are inverse of each other. f{x) appears in the 
equation of the trajectory in the {C,, () coordinates in the following way, 



e = /(C) 



(B.8) 



It is straightforward to check that the functions e*'"^*^"''"-' and e™^*-"'^^ are the two 
independent solutions of the Klein-Gordon fleld equation in the given metric. In 
our case the functions f{x) and g{x) appearing in eqs.(B.7) turn out to be f\{x) 
and g\{x) introduced in eqs.(4.32) and (A. 2) respectively so that the modes take 
the final form as shown in eqs.(4.45). 
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C Point-splitting Regular izat ion in the Draining Fermi 
Sea Background 

It will be useful for our purpose to work in the {u, v) coordinate system. The metric in 
this system is given by, 






where a,b = u, v. We get, 



T = 8Ae^9,09„0 + -^9„09,0 + -^9,09„0 . (CIO) 

Here we are distinguishing (9„0(9^0 from d^(j)du<p as we shall consider a point-splitting 
regularization of the operators. The components of the stress-tensor in the (m, v) system 
turn out to be, 

Tuu = <9„0(9„0 , 

Ty^ = d^(t)d^(t) + 2hx{vfdu<pdu(l) + hx(v){du<pd^(l) + d^(t)du<p), 

1 
Tuv = 2 (^u<pd^<P - d^(l)du4>) - hx{v)du(pdu(p , 

1 
Tyu = --{du(pdy(f) - d^(f)du(p) - hx{v)du(pdu(p , (Cll) 

where hx{v) has been defined in eqs.(5.55). At early time {v -^ — oo =^ hx{v) -^ 0) we 
recover the standard results for flat space. The vacuum expectation values can be written 
in terms of the basic objects, 

Dab{u,v) = {daM<P) . (C.12) 

For explicit computation we expand in terms of the in-modes. Then we point-split 
Dab{u,v) in the following way, 

D,t,{u,v) = lim / dw da,(t^Z'i^i,Vi)d,,<PZ*{u2,V2) , (C.13) 

OU,OV~fO JQ 

where the point-splitting is taken to be symmetrical (5.51). Using the full in-modes given 
in eqs.(4.45) and the result J^ dw we*'"^ = — ^ one gets. 



" ^A(t'i)/!;(t'2) ^ hx{v2)fx{vi) 



hxivi)hxiv2) f'xiVl)fxiv2) 



ifxiv2)-^iy ifxivi)-^2y 
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h\{v2) 

h\{vi) 
ATi{6if 



s 
A-n 

s 
An 



f'xiv2) 

(6 - fx{v2)r 

fxivi) 
(/a(^i)-6)^ 



(C.14) 



Notice that the terms in the square brackets do not have any short distance divergence. 
Moreover, because of the factor in the denominator, they are suppressed both at late time 
and at a finite time but large |g|. Therefore we shall simply drop these terms. The rest 
of the terms are computed, using the identity (5.61), to be. 



Dab — Sab + Fab , 

where the singular and finite terms Sab and Fab respectively are given by, 

1 






47r(7 - hxYiSv)^ ' 

hi 
47r(7 - hxY{5vY 
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hy 
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(C.15) 



(C.16) 
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27r(7 - hx)^ 167r(7 - hx^ 
hxi){hx) 



2A-K 



K 



27r(7 - hxf 327r(7 - hx)^ 



(C.17) 



where ip{hx) is given in eqs.(5.55). hx appearing on the right hand sides of all the above 
equations is meant to be evaluated at v. We have also used eq.(5.54). These results 
and the standard transformation laws of stress-tensor components under the coordinate 
transformation {u,v) -^ {t,q) allow us to finally compute {T^^). The results are, 

(T^.) = (T^.)5+(T^.)f, (C.18) 

where the singular and finite parts {T^u)s and {T^u)f respectively are given by. 
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(C.19) 
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S{fx,v) 
247r 



(C.20) 



Notice that each of {Tqi)s and (Tio)s has a hnearly divergent term which comes from the 
similar terms in S^v and S^u given in eqs.(C.16). 

D Relation between Bogolubov and Stress- Tensor Anal- 
ysis 

Some part of {T^y)^ can, in fact, be written in terms of the Bogolubov coefficients a{w, w') 
and l3{w,w'). Following [44] we expand the right hand side of eq.(C.12) in terms of the 
out-modes and use the unitarity properties of the Bogolubov coefficients to get. 



Dab 
Bab 



D:f + B, 



B 



(p) 

ab 



B. 



ab 1 

(m) 

ab 1 



(D.21) 



where 
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(D.22) 



Proceeding in the same way as the one that leads to eqs.(C.15) one can compute -0^^*. 
Subtracting this off from Dab given in eqs.(C.15) one gets Bab- The answers that one gets 
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are the following, 



1 1 

Bun = 7T7zSi9x,0 



247r '^ '^' 487r ' 

z47r 4o7r 

where we have taken the large time limit (lim„^oo S{g\,^) -^ —1/2) in the above expres- 
sions. Notice that Bab can also be computed using eqs.(D.21) and (D.22) in terms of the 
Bogolubov coefficients. This requires some particular integrals of these coefficients to be 
equal to functions appearing on the right hand side of eq.(D.23). It can be explicitly 
checked that all such relations can be satisfied by the following single equality. 

li'\v) + li^\v) = -^Sih,v), (D.24) 



where, 



Ix\v) = r dwidw2^M^e'^'"'~'"'>px{wi,W2) 
Jo 



l[^\v) = Re 






(D.25) 



For genuine Bogolubov coefficients satisfying the unitarity conditions the above relation 
must be true. Our Bogolubov coefficients given in eqs.(4.36) are approximately correct 
for small A. Although we have not quite succeeded, we attempted to understand how far 
the equality (D.24) can be satisfied by our Bogolubov coefficients. Below we shall report 
on how far we have been able to go. 

Using the variable y (5.70) one can erase the appearance of A on the right hand side 
of (D.24). 

S{f,,v) = -dyh{y)-^h{yf . (D.26) 

This means that it should be possible to fully absorb the A and v dependence of (/}'' (f ) + 
/|^ (f )) into y. Using the a and /3 coefficients in (4.36) one can readily check that it really 
happens for /} (f), 

I^P'iy) := r{'{v) = — / dwidw2 ^rr- — TFT" — V / ^^ ^T ' 

Att Jo 1 [tWiJL {—tW2) Jo smhTTW 

(D.27) 
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but not quite for I^^'{v), 

Viitiu 1 I- 



li^\v) = -i-Re 

Re 

An 



r{iwi 

rfiwijrfl + iw2) sinnTTW 



F(77/Ii iF M -I- 7'7/Io 1 tlin n 777/! ^. J 



(D.28) 



Notices that because of the factor kept in the curly brackets in the above expression one 
can set, for small A, 

li'\v) = 0, (D.29) 

so that one is left with the identity, 

I^'\y) = ^dyhiy) + ^Hyf . (D.30) 

It has not been possible for us to check this identity either analytically or numerically. 
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